What buoyancy really is. A Generalized Archimedes Principle 
for sedimentation and ultracentrifugation 
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Particle settling is a pervasive process in nature, and centrifugation is a much versatile separation technique. 
Yet, the results of settling and ultracentrifugation experiments often appear to contradict the very law on which 
they are based; Archimedes Principle - arguably, the oldest Physical Law. The purpose of this paper is delving 
at the very roots of the concept of buoyancy by means of a combined experimental-theoretical study on sedi- 
mentation profiles in colloidal mixtures. Our analysis shows that the standard Archimedes' principle is only a 
limiting approximation, valid for mesoscopic particles settling in a molecular fluid, and we provide a general 
expression for the actual buoyancy force. This "Generalized Archimedes Principle" accounts for unexpected 
effects, such as denser particles floating on top of a lighter fluid, which in fact we observe in our experiments. 



Sedimentation of particulate matter is ubiquitous in the nat- 
ural environment and widespread in industrial processes. For 
instance, particle and biomass settlingis responsible for the 
formation of depositional landforms yj] and plays a crucial 
role in marine ecology [2], while centrifugation of insolu- 
ble solids is a valuable separation methods in the extractive, 
chemical, and food processing industry 0]. Thanks to the ge- 
nius of Jean Perrin, sedimentation studies also provided the 
key support to the theory of Brownian motion [4], and origi- 
nated powerful methods to investigate soft and biological mat- 
ter, such as ultracentrifugation, a standard tool to obtain the 
size distribution of biological macromolecules or to pellet cel- 
lular organelles and viruses IHl. A particle settling in a simple 
fluid is subjected, besides to its weight, to an upward buoy- 
ancy force that, according to Archimedes' principle, is given 
by the weight of the displaced fluid. Usually, however, the 
settling process involves several dispersed species, either be- 
cause natural and industrial colloids display a large size distri- 
bution, or because additives are put in on purpose. The latter is 
the case of density-gradient ultracentrifugation (DGU), where 
heavy salts, compounds like iodixanol, or more recently col- 
loidal nanoparticles, are added to create a density gradient in 
the solvent. In DGU, proteins, nucleic acids, or cellular or- 
ganelles are expected to accumulate in a thin band around the 
position in the cell where the local solvent density matches 
the density of the fractionated species, the so-called isopycnic 
point. 

DGU is extremely sensitive, allowing for instance to re- 
solve differently labeled genomes with high efficiency |01, 
yet a subtle puzzle recurs in several studies. Even in earlier 
DGU measurements, the apparent density of some proteins 
was found to depend on the medium used to establish the 
density gradient [7]. The advent of sol-based DGU, allowing 
not only for more efficient separation of cells 1*81 IH, but also 
for fractionation of carbon nanotubes [10] and graphene [{Til, 
brought out more striking discrepancies. Indeed, the isopycnic 
densities of organelles |8] or carbon nanotubes fl^ fraction- 
ated using Percoll^*^, a standard DGU sol, are markedly dif- 
ferent from those found in sucrose or salt gradients, and strik- 
ing anomalies have been observed even for simple polystyrene 



latex particles Ill2ll . What value should then we take for the 
density of the medium, to predict the isopycnic point, if the 
surrounding fluid is not a simple liquid, but rather a com- 
plex mixture including other particulate species of different 
size and/or density? Similar ambiguities exist in experimental 
and numerical studies of colloid mixture settling in fluidized 
beds iflsilTill . where it is highly debated whether the density p 
of the bare solvent, or rather the density ps of the suspension 
should be used to evaluate the buoyant force. The latter choice 
is more widespread, but both attitudes have been taken in the 
literature [15], and even empirical interpola ting ex pressions 
have been suggested to fit experimental data IH 



The key point of our argumentation is that, when the sus- 
pending fluid is a colloidal suspension or a highly structured 
solvent, the amount of "displaced fluid" occuiTing in the sim- 
ple Archimedes' expression is substantially modified by the 
density perturbation induced by the particle itself in the sur- 
rounding. We shall focus on binary mixtures of particles of 
type 1 and 2, whose volumes and material densities are respec- 
tively given by (Vi, pi) and {V2, P2), suspended in a solvent 
of density p, under the assumption that component 1 is very 
diluted. Let us consider, as in Fig. [T] a large spherical cavity 
of volume V surrounding a single type-1 paiticle, and try to 
extend the common argument used to derive the Archimedes' 
principle. In the absence of paiticle 1, mechanical equilibrium 
requires the total pressure force exerted by the external fluid 
on V to balance exactly the weight W — m2n2gV, where 712 
is the number density of type-2 particles and 1112 — {p2~p)V2 
their buoyant mass. When particle 1 is inserted, however, the 
distribution of type-2 particles in V changes, because interac- 
tions generate a concentration profile set by the mutual radial 
distribution function 512 (r), which quantifies the local devia- 
tions from uniform density [18]. The total weight of the type- 
2 particles in V is now given by W' — m2gn2 /y gi2{r)d'^r. 
By taking the size of the cavity much larger than the range of 
512 (r), the total mass contained in V will then be subjected to 
an unbalanced mechanical force ll25ll 



Fi = W-W' = -magna / [512 (r) - 1] dV. (1) 
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FIG. 1: (Color on line) Schematic view of the density perturbation 
induced in the surrounding fluid by a settling colloidal particle or ra- 
dius Ri , with the upper panel showing the mutual radial correlation 
function gi2{r) of type-2 particles or radius R2 = qRi. The small 
q (or low-density) approximation leading to Eq. (|4) corresponds to 
evaluate density changes by just taking into account the white "de- 
pleted" spherical shell lying between Ri and -R2. 



Provided that the density correlations embodied by guir) are 
fully established, Fi will also amount to an effective excess 
buoyancy force acting on the test particle, which adds up to 
the usual Archimedes' term Fq = —pVig. This "Generalized 
Archimedes Principle" (GAP), which is our main theoretical 
result, can be equivalently written in terms of purely thermo- 
dynamic quantities. Provided that the number density rii of 
type-1 particles is very low, it is indeed easy to show that (see 
Supplementary Material): 
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where 11 is the osmotic pressure of the suspension. Eq. (|2) 
shows that Fi is proportional to the buoyant mass of type-2 
particles and to the osmotic compressibility, whereas the last 
factor explicitly accounts for mutual interactions between the 
two components. 

For spherical particles of radii Ri and R2, a simple expres- 
sion for Fi can be derived provided that component 2 is very 



diluted too, or, alternatively, that the range of gi2{r) is much 
smaller than which is usually the case if the size ratio q — 
R2/R1 < 1. In this limit, taking .912 (r) = Oforr < R1+R2, 
and 1 otherwise, we get Fi = (47r/3)(i?i + i?2)'^fi2"^2.'7- This 
result has a simple physical explanation: the excess buoyancy 
comes from the type-2 particle excluded from the depletion 
region shown in white in Fig.l. The total buoyancy Fi + Fq 
yields an "effective" density of the suspending fluid 



p* = P + $2(i + g)^(P2 -p), 



(3) 



where $2 is the volume fraction of type-2 particles. Note that, 
assuming p2 > p, p* is always larger than both p and ps — 
p+ {p2 — p)^2- Hence, the empirical interpolating expression 
suggested in lfl6ll is incorrect. A straightforward consequence 
is that the weight of a type-1 particle is exactly balanced by a 
suspension of type-2 particles at volume fraction: 



(4) 



which can be substantially lower than the isopycnic value 
'''2''° — (Pi ~ P)/(P2 — p) one would get from assuming p* 
equal to the suspension density. In the general, however, the 
additional force Fi may not necessarily oppose gravity. A 
strong attractive contribution to the mutual interaction may 
indeed overbalance the excluded volume term we considered, 
reversing the sign of Fi. Hence, particle 1 can actually be 
pulled down by the surrounding, showing an apparently larger 
density. 

Although derived for colloid mixtures, Eq. ([T) is valid in 
much wider conditions, whenever the region of perturbed sol- 
vent density is not negligible compared to Vi . Moreover, be- 
ing solely based on a force balance argument, Eq. ([T) does not 
require the suspension to have reached sedimentation equilib- 
rium, but only that the density distribution of type-2 particles 
around particle 1 has fully settled. Hence, since the time scale 
for the latter is usually much faster (at least for Brownian par- 
ticles), these predictions could be in principle checked on set- 
tling mixtures or in fluidized bed experiments. In practice, 
however, telling apart buoyancy effects from viscous forces is 
quite hard, because of the presence of long-range hydrody- 
namic interactions lll9ll . 

Thus, to test these ideas, we have devised a targeted equi- 
librium measurement. We have studied model colloidal mix- 
tures, obtained by adding a minute quantity ($1 < 10~^) 
of polymethyl-methacrylate (PMMA, pi — 1.19 g/cm'', ob- 
tained from microParticles GmbH, Berlin) particles with three 
different particle sizes (i?i ~ 220, 300, 400 nm), to a moder- 
ately concentrated suspension of spherical particles with ra- 
dius i?2 = 90 nm made of MFA, a tetrafluoroethylene copoly- 
mer with density p2 = 2.14g/cm^ ll20ll . MFA particles, 
though spherical and monodisperse, are partially crystalline, 
and therefore birefringent. Their intrinsic optical anisotropy 
yields a depolarized component Ivh in the scattered light 
that does not depend on interparticle interactions, but only 
on the local particle concentration [|20|] . Hence, the full equi- 
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FIG. 2: (Color on line) Inset: Equilibrium sedimentation profile of a 
suspension of MFA particles with radius R2 — 90 nm, dispersed in 
a solution of urea in water with density p = 1.04 g/cm'^. Here z is 
the distance from the cell bottom, $2 (z) the local MFA volume frac- 
tion, and the full line is the theoretical profile for hard-spheres with 
a radius R' ~ l.li?2. On the profile, the mean position of the thin 
layers of PMMA particles with radius 400 (bullet), 300 (triangle), 
and 220 nm (square) are compared to the prediction from the sim- 
ple Archimedes' principle (open dot, corresponding to $2 ~ 0.136). 
Main body: Expanded view of the profile region within the rectangu- 
lar box in the inset, showing the local density ps of the MFA suspen- 
sion. Superimposed are the full distributions (with normalized area) 
of the PMMA particles obtained from turbidity measurements and 
fitted with gaussian distributions as described in the text. Note the 
location of the isopycnic point where pa = p2- 



librium sedimentation profile can be simply obtained by ver- 
tically scanning a mildly focused laser beam and measuring 
IvH as a function of the distance from the cell bottom. A 
simple numerical integration of the experimental profile yields 
moreover the full equation of state of the system 112 U |22 
addition, MFA has a very low refractive index n 



In 

1.352, so 

it scatters very weakly in aqueous solvents. For better index- 
matching, we have used as solvent a solution of urea in water 
at 15% by weight, with density p — 1.04 g/cm'^. Hence, at 
equilibrium, the PMMA particles can be visually spotted as a 
thin whitish layer lying within a clear MFA sediment. 

The equilibrium sedimentation profile of the MFA sus- 
pension obtained by DeLS, is shown in the inset of Fig. |2l 
Using the simple Archimedes' principle, we would expect 
the PMMA particles to gather around the isopycnic level, 
namely, the region where the local suspension density is about 
1.19 g/cm'^, which corresponds to = 0.136. However, 
the layers lie well above this level, the more the smaller the 
PMMA particles are. The distribution of the guest parti- 
cles can be obtained by evaluating via turbidity measurements 
the sample extinction coefficient through the layer, where the 



TABLE I: Theoretical and experimental values for the effective 
isopycnic points $2 and for the standard deviation of the gaussian 
fits to the PMMA profiles. Calculated values are based on the simple 
"excluded volume" approximation leading to Eq. (|4j and (|5}, which 
may be reasonably expected to hold because the values of $2 are 
rather small and q not too large. 
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PMMA peak concentration does not exceed <I>i ~ 10~ . The 
body of Fig. |2]shows that the normalized probability distribu- 
tions for the PMMA particle position have a bell shape cen- 
tered on anomalously high z-values, with a width that grows 
with decreasing PMMA particle size. Since the MFA profile 
changes very smoothly on the scale of the layer thickness, it 
is in fact easy to show (see Supplementary Material) that the 
PMMA particles should approximately distribute as a gaus- 
sian with standard deviation: 
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(5) 



where tgi = ksT/ mig is the gravitational length of the type- 
1 particles, which we assume to be much larger than Ri and 
i?2. Table 1 shows that the experimental values agree very 
well with the values predicted by Eq. both for the effec- 
tive isopycnic point $2 and for the standard deviations of the 
gaussian fits. 

When considering the opposite case of small, dense par- 
ticles settling in a "sea" of larger but lighter ones, the GAP 
yields rather surprising predictions. Eq. (|2) shows indeed that 
Fi is proportional to the weight of a large particle: actually, 
the density perturbations in the host suspension can generate 
an excess buoyant force Fi amounting to a sizable fraction of 
mig, thus yielding an upward push on the small particle that 
largely outbalances its own weight. More specifically, in the 
Supplementary Material we show that, for hard-sphere mix- 
tures with q ^ 1, Fi is strongly non-monotonic, reaching 
a maximum at $2 < 0.2. Hence, most of the denser par- 
ticles will accumulate atop the lighter ones ll26ll . A striking 
example of this rather weird effect is shown in Fig. |3] where 
gold particles, with a radius of about 16 nm and a density 
pi ~ 19.3 g/cm'^ are seen to float mostly in the upper, very 
dilute region of an equilibrium sedimentation of MFA parti- 
cles (here q ~ 5.6). The DeLS profile shows that the MFA 
suspension is actually a colloidal ^M/ii (not a solid), with a 
density as low as ps ~ 1.2g/cm'' around the region where 
most of the gold particles accumulate. Since, for $2 ^ 0, 
the excess buoyant force Fi vanishes, some of the latter must 
lie within the MFA fluid phase too with a concentration pro- 
file that decreases downwards, as confirmed by turbidity data. 
Similarly, gold particles are expected to distribute in the su- 
pernatant solvent too, according to a barometric law c{z) cx 
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FIG. 3: (Color on line) Equilibrium sedimentation profile (A) and 
visual appearance (B) of a MFA suspension with a little amount of 
_Ri ~ 16 nm gold particles added. As evidenced by the weak Bragg 
reflections, the phase closer to the cell bottom is a colloidal crystal, 
whereas the upper phase is a colloidal fluid. To enhance the visibility 
of the thin gold layer, the sample picture has been taken using green 
illumination, with a narrow wave-length band around the plasmonic 
absorbtion peak of the gold colloid. The concentration profile ob- 
tained from turbidity data (exploiting in this case the proportionality 
between gold absorption and local concentration) shows that gold 
particles are also present both within the MFA sediment and in the 
supernatant solvent. The semilog plot of the polarized scattering in- 
tensity in Panel C is fitted with a double exponential, as discussed in 
the text. 



exp(—z/£gi), with a gravitational length £gi ~ 1.4mm. This 
weak barometric region can be detected by polarized light 
scattering 1,271 . Panel C in Fig. |3] shows that the polarized 
scattered intensity can be fitted as the weighted sum of two ex- 
ponentials I ~ Ii exp{—z /£gi) + 12 exp{~z/£g2), where the 
MFA gravitational length is fixed at the value £g2 — 0.13 mm, 
whereas from the fit £gi ~ 1.38 mm for gold. This value for 
£gi corresponds to an average particle radius Ri ~ 16 nm, in 
very good agreement with the estimate made from the position 
of the particle plasmonic absorption peak at A = 528 nm. 

The GAP qualitatively accounts for the anomalous DGU 
measurements of polystyrene bead density [12J, even when, 
in the presence of oppositely-charged nanoparticles, the lat- 
ter apparently increases, and for em piri cal expressions used 
to fit flotation-bed experiments ilH IitIi. But Eq. ([T]i has a 
much wider scope. For instance, provided that a model for 
gi2 is available, it should correctly account for "solvation" ef- 
fects on the buoyancy force felt by proteins, simple molecules, 
even single ions, or provide a sensitive way to detect by DGU 
aggregation and association effects in biological fluids. Sim- 
ilarly, corrections to the simple Archimedes' expression will 
also show up for nanoparticles settling in a strongly corre- 
lated solvent, such as a pure fluid or a liquid mixture close 
to a critical point. Some relation with the Brazil nut effect 



in granular fluids, which is also affected by the densities of 
the grain |23, 24], may also exist, although the latter is usu- 
ally complicated by the presence of dissipation, convective ef- 
fects, and effective thermal inhomogeneity. In fact, due to 
its exquisite sensitivity to the specific properties of a mixture, 
the "reversed" gravity-segregation effect we have highlighted 
may allow to devise novel sophisticated DGU fractionation 
methods, able to tell apart solutes with the same density and 
composition, but different size. 
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[25] We assume that the both mutual interactions between the two 
species and self interactions between type-2 particles are suf- 
ficiently short-ranged. Eq. ([T) is also valid when the host sus- 
pension is non uniform, provided that 712 varies slowly over the 
range of 312(7-). 

[26] Note that the accumulation on top of the heavier particles does 
not however lead to a macroscopically inverted density pro- 
file. An attentive examination of Eq. l[T) shows indeed that the 
weight increase with respect to a suspension of type-2 parti- 
cles at volume fraction <I>*, due to the presence of the heavier 
particles, is exactly balanced by the "expulsion" from the accu- 
mulation layer of those particles of type 2 that yield the excess 
buoyancy F\ . Macroscopic hydrodynamic stability is thus pre- 
served. 

[27] Although in index-matching, MFA particles still scatter polar- 
ized light, which is however fully incoherent and proportional 
to the depolarized scattered intensity ||20|| . 



compressibility sum rules (Ref. [17]): 



n2 S{G) =kBT 
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According to the McMillan-Mayer theory of solutions, the contribu- 
tion of the solvent to the total free energy can be disregarded if effec- 
tive interactions among particles are introduced. In the limit 711 — )■ 
we can express the free energy derivatives appearing in Eq. (|S5t in 
terms of the osmotic pressure: 
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leading to 



SUPPLEMENTARY MATERIAL 

Effective buoyancy. 

We provide here a formal derivation of the buoyancy force Fi act- 
ing onto a test type-1 colloid immersed in a solution of type-2 parti- 
cles, expressed in purely thermodynamic terms. The density profile 
of a suspension of particles in the presence of a gravitational field is 
described by the hydrostatic equilibrium condition 



dU[n2{z),T] 
dz 



-m2gn2[z) 



(SI) 



where m2 is the buoyant mass of type-2 particles, 11 the osmotic 
pressure, and we assume that the number density n2 may depend on 
z. The gravitational length £g = ksT/ (m2g) defines the charac- 
teristic scale of the spatial modulations of the density profile: here 
and in the following we assume that £g is the largest length in the 
problem, a condition easily met in colloidal suspensions. Under this 
assumption, the contribution to the buoyancy force acting onto a test 
particle (denoted by index 1) inserted in this solution, due to the pres- 
ence of type-2 particles, is given by Eq. (1) in the text: 



fi (2) = —7712 (? 712(2) / drhi2ir) 
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(S2) 



where hi2{r) — guij) — 1. This expression depends on the mutual 
correlations between the two species but can be equivalently written 
in terms of purely thermodynamic quantities. Regarding the sys- 
tem as a binary mixture where component 1 is extremely diluted, the 
Ornstein-Zemike relation in the n\ Q limit (see Ref. [14]) 

/ii2(7-) = Ci2(r) + 712 y dxci2(r - x) /l22(x) (S3) 

allows to express the integral of h\2 {r ) in terms of the integral of the 
direct correlation function C12 (r) and the long wave-length limit of 
the structure factor of a type-2 one component fluid 5*22(0): 
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which coincides with Eq. (2) in the paper. This shows that the contri- 
bution to the buoyancy force on a type-1 particle due to the presence 
of component 2 is proportional to the buoyant mass 7712. It is interest- 
ing to investigate the limiting form of the buoyancy force when the 
type-1 particle is just a "tagged" type-2 particle, with identical phys- 
ical properties. In this case the system is effectively one-component 
and then = . The buoyancy force acting onto a particle in 
the solution acquires the form: 



F — rag 



l-kgT 



m 

dn 



(S8) 



It is instructive to deduce Eq. ( IS8t with a different approach, 
which highlights its physical meaning. The equilibrium sedimenta- 
tion profile of a suspension of interacting Brownian particles is usu- 
ally derived by balancing gravity with the diffusive term deriving 
from gradients in the osmotic pressure. However, fixing the attention 
on a single test particle, we can try to summarize the effect of all the 
other particles as an "effective field" F adding to the bare gravita- 
tional force —mg. From the Smoluchowski equation, the combina- 
tion of these two contributions yield a density profile: 

, „dn , ^ 
ksT— ^n{F -mg), 
dz 

that, combined with the hydrostatic equilibrium equation JSlt . yields 
for F the expression in Eq. ( |S8) . Hence, the equilibrium sedimenta- 
tion profile of an interacting suspension can be equivalently viewed 
in terms of the probability distribution for the position of a test parti- 
cle subjected to a spatially-varying gravitational field, whose depen- 
dence on z is dictated by the equation of state of the suspension. 

In hard sphere systems we can easily obtain an approximate ex- 
pression for the buoyancy force from Eq. ( IS71 : a rough estimate of 
the excluded volume effects in the osmotic pressure can be obtained 
following the familiar Van der Waals argument: 



drhi2{r) = S22(0) / drci2(r-) 



(S4) 



Both terms at right hand side can be expressed as thermodynamic 
derivatives of the Helmholtz free energy of the mixture A via the 
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By substituting this form into Eq. l lS7t we recover tlie simple result, 
already quoted in a slightly different form in the main paper 

i^i = m2g<S>2 ^1 + (SIO) 

A more careful evaluation is obtained by starting from the analytical 
expression of the excess free energy of a binary hard sphere mixtures 



provided by Mansoori et al. (J. Chem. Phys. 54, 1523, 1971). The 
result can be conveniently expressed in terms of the effective mass 
density of the surrounding medium p* defined by 

Fi^^TvRlp*g (Sll) 
The explicit expression for the effective density reads: 



p* ^ 6 + (1 - qf{2 + g) (1 - $2)^ - 3(1 - g^)(l - $2)" - 2 [(1 - qf{2 + q) - q'] (1 ~ <i>2) 
m2n2 (1 - $2)-* + $2(8 - 2$2) 

I 



The dependence of p* on the size and volume fraction of type-2 par- 
ticles is shown in Fig. |4] For g > 1, i.e. when a small test particle 
is immersed into a suspension of big particles, the buoyancy force 
displays a pronounced maximum. In the g — >^ 00 limit, the maxi- 
mum buoyancy force is attained at <I>2 ~ 0.154, where it reduces 
to a sizeable fraction of the effective weight of a type-2 particle: 
Fi ~ 0.055 maff. 




FIG. 4: Effective mass density of the surrounding medium, relative to 
the type-2 mass density as a function of $2 for different g = R2/R1. 
Left panel: results for g < 1. Right panel: g > 1. Note the change 
of scale in the vertical axis. 



Distribution of guest particles at equilibrium. 

The hydrostatic equilibrium condition for a suspension of type-1 
particles reads: 

^=ni [-mig + Fi] (SI2) 
az 

where 11, ni and mi are the osmotic pressure, average local den- 
sity and buoyant mass respectively. In the limit of short range inter- 
species correlations, the excess buoyant force Fi due to the presence 
of type-2 particles is given by Eq. ( ISIOt , while in the diluted limit of 
type-1 particles the ideal gas equation of state Hi = nifcsT holds. 
Substituting these results in Eq. ( ISI2) we find: 



ks-i — — =nig 
dz 



-nil + 1712 $2(2:) 1 + 



1 



(S13) 



which defines the number density profile of type-1 particles. The 
maximum of the resulting distribution corresponds to the vanishing 
of the right hand side of this expression, given by condition (3) of the 
main paper: 



$2(2*) = 



(S14) 



{m\/m2)q^ coincides with the isopycnic volume 
fraction defined in the main paper. By expanding $2 (z) around the 
position of this maximum z*, Eq. l |SI3t becomes: 



dn\ m2 q d^2(z) 

=ni 

dz ksT dz 



d$2{z) 
dz 



e,i'i>2 



(SI4) 



whose solution ni (z) is a gaussian centered in z = z* with standard 
deviation given by Eq. (5) of the main paper. 



